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1. Introduction 

Approximation of set-valued functions (SVFs) has various potential applications in optimization, control 
theory, mathematical economics, medical imaging and more. The problem is closely related to the approxi- 
mation of an iV-dimensional object from a sequence of its parallel cross-sections, since such an object can be 
regarded as a univariate set- valued function with sets of dimension A'^ — 1 as images [11, 25]. In particular 
for = 3, the problem is important in medical imaging and is known as "reconstruction from parallel 
cross-sections" (see e.g.[l, 5, 6] and references therein). 

Motivated by the problem of set-valued approximation and its applications, we develop and study set- 
valued subdivision schemes. Real-valued subdivision schemes repeatedly refine numbers and generate limit 
functions. When applied componentwise to points in ]S?, the schemes generate smooth curves/surfaces, and 
as such, are widely used in Computer Graphics and Geometric Design. When the initial data are samples of a 
function, the limit of the subdivision approximates the sampled function. For a general review on subdivision 
schemes see [14]. In this work, we propose a new method for the adaptation of subdivision schemes to sets, 
and show convergence and approximation properties of the resulting set-valued subdivision schemes. 

In the case of data consisting of convex sets, methods based on the classical Minkowski sum of sets can 
be used [9, 30]. In this approach, sums of numbers in positive operators for real- valued approximation are 
replaced by Minkowski sums of sets. A more recent approach is to embed the given convex sets into the 
Banach space of directed sets [3] , and to apply any existing method for approximation in Banach spaces [4] . 

The case of data consisting of general sets (not necessarily convex), which is relevant in many applications, 
is more challenging. For data sampled from a set-valued function with general sets as images, methods based 
on Minkowski sum of sets fail to approximate the sampled function [12, 30]. So other operations between 
sets are needed. 

In the spirit of Frechet expectation [18], Z. Artstein proposed in [2] to interpolate data sampled from a 
set- valued function in a piecewise way, so that for two sets Fi,F2 C M" given at consecutive points a;i,a;2, 
the interpolant F {■) satisfies for any ^1,^2 G [2^ii2^2]i 

dist(F{ti),F(t2)) - ^^^~^^^ dist (fl,f2), (1.1) 
X2 - Xi 

where dist is a metric on sets. The relation (1.1) is termed the metric property. A weighted average of two 
sets introduced in [2] , and termed later as the metric average, leads to a piecewise interpolant satisfying the 
metric property relative to the Hausdorff metric. 

Extending this work, in [10] the metric average is applied in the Lane-Riesenfeld algorithm for spline 
subdivision schemes. The set-valued subdivision schemes obtained this way are shown to approximate SVFs 
with general sets as images. The convergence and the approximation results obtained in [10] are based on 
the metric property of the metric average. The adaptation to sets of certain positive linear operators based 
on the metric average is described in [16]. For reviews on set-valued approximation methods see also [13, 27]. 
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Fig 1. The graph of a SVF F (x), which is Lipschitz continuous relative to the symmetric difference metric, while discontinuous 
relative to the Hausdorff metric. 

As it is noticed in [2], the particular choice of metric is crucial to the construction and analysis of set- 
valued methods. While previous works develop and analyze set-valued approximation methods in the metric 
space of compact sets endowed with the Hausdorff metric, we consider the problem in the metric space of 
Lebesgue measurable sets with the symmetric difference metric^. Our setting allows us to approximate SVFs, 
which are Holder continuous in the symmetric difference metric but may be discontinuous in the Hausdorff 
metric, as illustrated by the following simple example. 

Example 1.1. Let F he a SVF from M to subsets o/M, 

F(x) = {2/:l<y<2-|x|}|J{y:3<j/<4-2|x|} , (1.2) 

with graph given in Figure 1. It is easy to observe that F is discontinues at x — (and also at x = —1,-^,1^, 
if the distance between subsets of M is measured in the Hausdorff metric, but it is Lipschitz continuous 
everywhere if the distance is measured in the symmetric difference metric. 

On the other hand, under mild assumptions on the sets F (x), Holder continuity in the Hausdorff metric 
implies Holder continuity in the symmetric difference metric [19]. Since there is an intrinsic connection 
between the continuity and approximability of a function, approximation results for Holder continuous SVFs 
obtained in the symmetric difference metric apply to a wider class of functions than similar results obtained 
in the Hausdorff metric. 

In order to develop approximation methods in the space of sets with the symmetric difference metric, we 
introduce a new binary weighted average of regular compact Jordan measurable sets. The new average is built 
upon the method introduced in [25], which is known as the shaped-based interpolation in the engineering 
literature [21, 28]. Our new average has the metric property relative to the symmetric difference metric. In 
addition when both weights are non-negative, the measure of the average of the two sets is equal to the 
average with the same weights of the measures of the two sets. We term this feature of the new average the 
measure property, and term our average the measure average. 

The measure average performs locally on each connected component of the symmetric difference of the 
two operand sets, leading to satisfactory geometric performance, which is essential in many applications. In 
particular, the ideas of this work lead to a practical algorithm for the reconstruction of 3D objects from their 
2D cross-sections described in [23]. 

First we use the measure average to interpolate between a sequence of sets in a piecewise way. Then we 
adapt to sets spline subdivision schemes, expressed in terms of repeated binary averages of numbers using 

^The measure of the symmetric difference is only a pseudo-metric on Lebesgue measurable sets. The metric space is obtained 
in a standard way as described in Section 2 
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the Lane-Riesenfeld algorithm [24]. As in the case of the metric average [11], we prove convergence of the 
sphne subdivision to a Lipschitz continuous hmit SVF F°° (•). It follows from the measure property of the 
measure average, that /i (•)) is the limit of the same spline subdivision scheme applied to the measures 
of the initial sets. Moreover, we prove that spline subdivision schemes adapted to sets with the measure 
average are monotonicity preserving in the sense of the set-inclusion relation. 

It is well known, that in order to obtain "reasonable" interpolation methods other than the piecewise 
interpolation, some notion of extrapolation is needed. Our measure average of sets is defined for positive 
weights and also when one weight is negative, therefore it performs both interpolation and extrapolation. 
Using the measure average with both negative and non-negative weights, we adapt to sets the 4-point 
interpolatory subdivision scheme of [15]. This is the first adaptation of an interpolatory subdivision scheme 
to sets. We prove that the 4-point subdivision scheme based on the measure average converges to a continuous 
limit SVF and approximates Holder continuous SVFs, when the initial sets are samples of such a function. 

We observe that many results on set-valued subdivision obtained in this and previous works are based 
on the triangle inequality in the underlying metric space, along with the metric property of the average of 
sets. Employing this observation, we extend several results obtained in the context of metric spaces of sets, 
to general metric spaces endowed with an average having the metric property. 

The structure of this work is as follows. Preliminary definitions are given in Section 2. In Section 3, we 
study properties of the average of sets based on the method in [25] , which are relevant to the construction 
of our measure average. In Section 4 we introduce our measure average of sets, prove its important features 
and apply it to the interpolation between sets in a piecewise way. Spline subdivision schemes based on 
the measure average are studied in Section 5, while in Section 6, we adapt to sets the 4-point subdivision 
scheme. In Section 7, we provide several computational examples illustrating our analytical results. Finally, 
extensions of some of the results obtained in metric spaces of sets to general metric spaces are discussed in 
Section 8. 

2. Preliminaries 

First we introduce some definitions and notation. We denote by fi the n-dimensional Lebesgue measure and 
by £„ the collection of Lebesgue measurable subsets of M" having finite measure. The set difference of two 
sets A, B is, 

A\B = {p:p(^A,piB} , 
and the symmetric difference is defined by, 

AAB = A\b|Jb\A . 

The measure of the symmetric difference of A, S e £„, 

d^ {A, B)^fi (AAB) , 

induces a pseudo-metric on £„, and d^) is a complete metric space by regarding any two sets A, B such 
that ^ {AAB) = as equal ([20], Chapter 8). For A, J5 e £„, such that B C A, it is easy to observe that 

d^{A,B)^piA\B) = p{A)-p{B) . (2.1) 

The boundary of a set A is denoted by dA, and we use the notation ci (A) for the closure of the interior of 
A. A bounded set A, such that A = ci (A) is called regular compact. Regular compact sets are closed under 
finite unions, but not under finite intersections, yet for A, B regular compact such that B d A, it holds 
trivially that An B = B = ci{Ar\ B). 

We recall that a set A is Jordan measurable if and only if p (dA) = 0. It is easy to see that for a Jordan 
measurable A, 

^l{A)^^i{ci{A)) . (2.2) 

We denote by 3n the subset of £„ consisting of regular compact Jordan measurable sets. Notice that for 
any A, B (z ^n, {A, B) ^ implies A ^ B, therefore d^ is a metric of J„. Moreover, by its definition Zn is 
closed under finite unions. 
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We recall that a set A is called connected if there are no two disjoint open sets Vi,V2 C M", such that 
A — n ^i) U n ^2)- The set C C D is called a connected component of D if it is connected, and if there 
is no connected set B, such that C C B C D. 

3. The "distance average" of sets 

The basic tool for the construction of the measure average of sets to be introduced in the next section is 
what we call the distance average of sets. In this section we derive properties of the distance average that 
are relevant to our construction. 

3.1. Definition and basic properties 

The distance average is based on the method introduced in [25] , which employs the signed distance functions 
of sets. The signed distance from a point p to a non-empty set A C M" is defined by, 

r d{p,dA) peA 

''^'^^^^\~dip,dA) Pi A, ^'-'^ 

where d {q, B) is the Euclidean distance from a point 5 to a set i?, namely 

d [q, B) — Ta\u\\q — b\\ . 
beB 

The signed distance function of A is defined on M" as dg (•, A). 

Definition 3.1. The distance average with the averaging parameter a; G M 0/ two not-empty sets A, B 
is, 

xA^il-x)B = {p: fA^BAp)>0} , (3.2) 

where 

fA,B.x (p) = xds (p, A) + il- x)ds [p, B) . (3.3) 

Note that fA.B,x is not the signed distance function of the set xA^ (1 — x) B. Also note that fA,B.x is 
continuous by the continuity of the distance function. 

We observe a few properties of the distance average that are relevant to the construction of the measure 
average in the next section. 

Lemma 3.2. Let A,B(^ 3:„ and x eM., then 

1. = B, 1A@QB = A 

2. xJ^{l-x)A^A ^ ^ 

3. For BCA,xi< X2, xiA^l - xi) B C X2A^ (1 - X2) B 

4. Forjre [0,1], AC\B C xA^{l- x)B C A\JB 

5. xA^ {1 — x) B is a bounded closed set 

Proof. Properties 1-2 follow from Definition 3.1. To obtain Property 3, observe that ds{p, B) < ds{p, A) and 
consequently fA,B,xt (p) < fA,B,x2 {p)- 

To prove Property 4, note that for p e A{^B, ds {p,A) > 0,ds ip,B) > 0, so for x e [0, 1], fA,B,x > 
and p e xAQ{l — x)B. Now with a similar argument for p ^ A[JB, Property 4 is proved. 

The set xA@ {\ ~ x) B \s closed by definition, so in order to prove Property 5, it is sufficient to show that 
for any x this set is bounded. Since A, B are bounded, so is their union. Therefore for x € [0, 1], Property 5 
follows from Property 4. For x ^ [0, 1], in view of Definition 3.1 we have, 



xA0 {l-x)B(^A[jB[j Ca.b.: 
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where 

Ca,b,x ^{p:ds (p, A) < 0, ds b, B) < 0, xds {p, A) + {I - x) ds {p, B) > 0} . (3.4) 

Since ^ IJ i? is bounded, it is enough to show that also Ca.b,x is bounded. Without loss of generality, assume 
that a; > 1. There exists > 0, such that for any p € M", 

\dsip,A)-dsip,B)\<9 . (3.5) 

From (3.4) and (3.5), 

Ca^b,. Q{p-ds {p, A) < 0, ds {p, B) < 0, xds {p, ^) + (1 - x) {ds {p, A) - 0) > 0} . (3.6) 
and therefore, 

Cam,x Q{p--ds {p, A) < 0, ds {p, B) < 0, Ids {P, A)\<{x-l)e} , 

from which it follows that Ca,b.x is bounded. □ □ 

We extend the domain of the distance average to include the empty set. In case B = cj), A cj), choose a 
"center" point q oi A such that, 

dsiq,A)^snp{ds{a,A)} . (3.7) 

aeA 

We define xA^{l — x)B as the set, 

{p:xds{p,A) + {x~l)\\p^q\\>0} . 

The average of two empty sets is the empty set. One can verify that with these definitions Properties 1-5 of 
Lemma 3.2 are preserved. 

The distance average does not satisfy the metric and the measure properties. For example, for any non- 
empty A,B G -3„ that are disjoint , ^A^^B is the empty set. Also consider the distance average ^A^^B 

with A ^ [0,3] and B = [0, 1]U[2,3]. From Property 4 of Lemma 3.2, we have that,i^0iB C A, on 
the other hand we now show the opposite inclusion. Let p G A \ B = (1,2), it is easy to observe that 
ds {p,A) > 0, ds {p,B) < and \ds {p,A)\ > \ds {p,B)\. Consequently, from the definition of the distance 
average p £ ^A^^B, thus ^A^^B = A. Of course this average does not satisfy the metric property or 
the measure property. Moreover, it is undesirable to obtain one of the original sets, as an equally weighted 
average of the two different sets, since such average does not reflect a gradual transition between the two 
sets. We aim to define a new set average with desirable properties, using the distance average as the basic 
tool for the construction. 

Before presenting the new set average, we consider the measure of the distance average as a function of 
the averaging parameter, 

h{x) ^ n(^xA^{l- x)B^ , X eR , (3.8) 

and study conditions for its continuity. Note that by Property 5 in Lemma 3.2, h (x) is well defined for all 

x€R. 

3.2. Continuity of the measure of the distance average 

First we prove a result, which might be of interest beyond its application in our context. 
Lemma 3.3. Let A, B he closed sets, and Zet A > 0, A 7^ 1. Define 

MA,B.x^[p:piA[jB,d{p,A) = Xd{p,B)] , (3.9) 

then n{Mji B.\) = 0. 
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Proof. To prove the claim of the lemma, it is sufficient to show that for any p £ Ma,b,\, there exists a 
cone of constant angle with p as its vertex, which is not in Ma,b.\- Without loss of generality assume that 
A > 1. Let p € Ma,b,\, then there exists an open ball of radius d(p,A) — \d{p,B) around p that contains 
no points of A, and another open ball around p, of radius d(p,B), which contains no points of B. There 
is at least one point v G B, such that — ?;|| = d{p,B). Let e e {0, d {p, B)) and be the point on the 
segment [p, u], such that ||a;e — v\\ — e. The open ball of radius Xd{p, B) — e around contains no points 
of A, and d {x^, B) ~ d (p, B) — e. Since A (d [p, B) — e) < \d {p, B) — e, there is no point of A at distance 
\d{x^,B) from Xj, and therefore x^ ^ Ma,b,x- Consider now a point x'^ at distance r from x^. By the 
triangle inequality, d{x'^,B) < d {x^, B) + r and d{x'^,A) > d {x^, A) — r. Therefore, 



d{x'^,B) <d{p,B)-e + r , (3.10) 



and 



d{x',,A)>Xd{p,B)~e-r. (3.11) 
In order to obtain that x'^ ^ Ma,b,\, it is enough to show that, 

Xd{x'^,B) <d{x'^,A) , (3.12) 

or using (3.10) and (3.11), 

Ad (p, B) - Xe + Xr < Xd{p, B) ~ e - r . (3.13) 

Solving (3.13) for r, we obtain. 

So the open ball of radius j^s around x^ has empty intersection with the set Ma,b,x- Let Up be the 
union of the balls B^ for e G (0, d {p, B)), 

Up^ y B,. (3.15) 

ee(0,d(p,S)) 

Any point q in the open ball of radius d{p, B) around p, such that the angle Zqpv satisfies, 

< tan Zqpv < ^— ^ , (3.16) 
A + 1 

is in Up. We conclude that for any p G M there is a cone Up of a constant angle based at p, such that 

U f] Ma.b.x = 0- In view of the Lebesgue density theorem {see. e.g. [26], Corollary 2.14), we obtain that 

Ma.b.x has zero Lebesgue measure. □ □ 

Note that it is easy to construct closed sets A, B with non-empty intersection, such that the claim of 
Lemma 3.3 does not hold for A = 1. 
From Lemma 3.3 we conclude, 

Corollary 3.4. Let A,B G and t eR. Define the set ^Ia.b.x o,s, 

^a^b.x = {p : fAM.x (p) = 0} . (3.17) 

Then for x ^, 



Proof. For x ^ 0,1, set A = 



(f^AB,^) = . (3.18) 

, then 

^A,B,x C MaA,dBA\JdA\JdB , 



with Mqamb.x defined by (3.9). It follows from Lemma 3.3 and the assumptions A, B G 3„ and x ^ ^, that 
fj,{ilA.B,x) = 0. For a: = or a; = 1, ^Ia.b.x equals to dB or dA respectively, so fJ.{^A,B.x) = 0, by the 
assumption that A,B E Zn- D n 
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Corollary 3.4, does not treat the case x = ^. Indeed, it is not difficult to give two sets A,B G such 
that /i ^Sl^ -B 5) Figure 2 for an example of such two sets. 

Since for the continuity of h{x) at x = ^, we need the condition fi(^flj^Q 1^ = 0, we introduce the 
following relation between two sets. 

Definition 3.5. A, B Cz 3n satisfy the zero-measure condition if ^ (^^a si) =0. 
Remark 3.6. In view of Corollary 3.4, if A, B satisfy the zero-measure condition, then 

(f^AB.x) = , (3.19) 

for any x. 

Remark 3.7. For a Lehesgue measurable function f : M" — > M, the set of values c such that, 

M({p:/(p)=c})>0, 

has zero measure in view of Fubini's theorem (see e.g. [20], Section 36). Consequently if the level set 
|-P ■ /a B i (-P) — o| ^'^■^ non-zero measure, one can always choose an arbitrarily small £ > such that, 

(^^P '■ Jab ^{p) ~ s|) = 0. Therefore, the case of A, B satisfying the zero-measure condition is generic, 
while the case that this condition is not satisfied is degenerate. 

Next we adapt to our context a basic result from probability theory. 

Lemma 3.8. Let A,B^ B C A and (JIa,b,x') — 0, then the function h{x) defined by (3.8) is 
continuous at x* . 

Proof. To sec that h{x) is left-continuous at x* , let Xn ^ x* , Xn < Xn+i. By Property 3 of Lemma 3.2, 

x„A0 (1 -Xn)BC a;„+iA0 (1 - Xn+i) B , 

so 

x*A (1 ~ x*) S - || Q XnA^ (1 - x„) B^ y nA,B,x' ■ 

Consequently by the continuity from below of the Lebesgue measure (see e.g. [22], Chapter 2, Section B, 
M5) and by the assumption fi {rtA.B,x') = 0, 

h (x*) = lim h {x„) . 

To obtain that h{x) is right-continuous at x* , let a;„ x, Xn > Xn^i. Then in view of Property 3 of Lemma 
3.2, 

x*A0 {l-x*)B=f] XnA (1 - x„) B . 

71=1 

Using Property 5 of Lemma 3.2, by the continuity from above of the Lebesgue measure {see e.g. [22], Chapter 
2, Section B, M6) we obtain that h{x) is right-continuous at x*. □ □ 

In view of Corollary 3.4, by Remark 3.6 and Lemma 3.8, we arrive at 

Corollary 3.9. If A, B G 3„ satisfy the zero-measure condition, then h(x) is continuous everywhere. Oth- 
erwise h{x) is continuous at all x except at x — |. 

Finally, we discuss the Jordan measurability of xA^{l — x)B, 

Corollary 3.10. If A, B G Zn satisfy the zero-measure condition, then xA^(l — x)B is Jordan measurable 
for all X. Otherwise it is Jordan measurable for all x except may be for x ~ ^. 
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Fig 2. The pentagon and the lower triangle represent the boundaries of the sets A and B respectively. The black set is contained 



Proof. By Definition 3.1, 

d [xA@ {x -t)B^ C nA,B,x , (3.20) 

which in view of Corollary 3.4 and Remark 3.6, leads to the claim of the corollary. □ □ 

However the distance average tA^ {1 — t) B is not necessarily regular compact, even when A and B are. 

4. Construction of the "measure average" of sets 

Our construction aims to achieve several important properties of the measure average, denoted hy tA ^ {1 — t) B . 
We use these properties later on in the analysis of the subdivision methods based on the measure average. 
For A, B £ and an averaging parameter t e M, the desired properties of the measure average are 

List of Properties 4.1. 

1. tAQ {1 — t) B £ (closure property) 

2. OA0 li? = B, lA^OB = A (interpolation property) 

3. df, {sA (1 - s) B, (1 - t) B) <\t- s\ {A, B) (suhmetric property) 

In addition for s,t Cz [0, 1], 

4. (tA0 (1 -t)B) ^ tfi (A) + (l-t) fj, (B) (measure property) 

5. d^ {sA (1 - s) B, (1 - t) B) = \t- s\ d^ {A, B) (metric property) 

6. ci (A fl B) C (1 - t) S C A U S (inclusion property) 

The above properties are analogous to those of weighted averages between non-negative numbers, defined 
by max{0, tp+ [1 — t) q], for p,q G M+ and i S M. In this analogy, the measure of a set, the measure of 
the symmetric difference of two sets (d^ {■,■)), and the relation C are replaced by the absolute value of a 
number, the absolute value of the difference of two numbers and the relation < respectively. 

The measure average is constructed in three steps, each based on the previous. 

4-1- The measure average of "simply different" sets 

We begin with the simple case oi A, B e such that B C A and A \ B consists of only one connected 
component. We call two such sets simply different. First we assume that the sets A, B satisfy the zero- 
measure condition. In this case, the measure average is a reparametrization of the distance average, so that 
the measure of the resulting set is as close as possible to the average of the measures of A, B, 

Definition 4.2. Let A, B be simply different sets satisfying the zero-measure condition. The measure average 
of A, B with the averaging parameter t e M is, 
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tA^{l-t)B = ci(g{t)A^{l-g{t))B] , (4.1) 
where g{t) is any parameter in the collection, 



X : X = a.rg mm \h{^) ~ {tfi {A) + {1 - t) iJ,{B))\\ . (4.2) 

i<£[-N,N] 

Here h is defined in (3.8) and N » 1 is a large positive number. 

The zero-measure condition satisfied by A, B and Corollary 3.9 imply the continuity of h. Therefore the 
argmin in (4.2) is well defined. By Corollary 3.10 and the relation (2.2), 

M(M0(l-t)i?) =/i(5(0) . (4.3) 

The measure average in Definition 4.2 is a regular compact set. Since h{xi) — h{x2), for any two parameters 
xi,X2 in the collection (4.2), the corresponding averages defined using either g(t) — xi or g{t) — xi in (4.1) 
have the same measure. In addition, by Property 3 of Lemma 3.2, one of these sets is necessarily contained 
in the other. Therefore, in view of (2.1) and because both sets are regular compact, 

a:iA0 (1 -x^)B^ 2:2^0 (1 -X2)B . 
So any parameter in the collection (4.2) can be used in (4.1). 

Lemma 4.3. The measure average of simply different sets satisfying the zero-measure condition has Prop- 
erties 1-6 in List of Properties 4^.1. 

Proof. The closure property follows from Definition 4.2. The interpolation property follows from Definition 
4.2 and from Property 1 in Lemma 3.2. 

Notice that since B <Z A^ Property 3 of Lemma 3.2 implies that h{x) defined in (3.8) is monotone non- 
decreasing. To obtain the submetric property, we denote m^_s = h{—N) and M^^b — h{N), where [—N, N] 
is the domain used in (4.2). Since the average ^ is a reparametrization of the average 

(1 - t) b) e [mA,B,MA.B] ■ (4.4) 



By the continuity of h, we obtain from Definition 4.2, that for any t satisfying, 

mA,B < tn {A) + {l-t)pi {B) < Ma,b , (4.5) 

we have, 

M {tA (1 - t) s) = t^i (A) + (1 - i) /i (B) . (4.6) 
Summarizing (4.4) and (4.6), we arrive at 

r ruA.B, tfi{A) + (l-t)fi{B)<mA^B 

//(i^0(l-t)s) - i t/i (A) + (1 - t) (S) , tfi(A) + (l-t)fi{B) e {mA,B,MA,B) (4.7) 

[ Ma,b, tfi(A) + (l-t)fi{B)>MA.B 

Assume without loss of generality that s < i, then from the monotonicity of h and Definition 4.2, g{s) < g(t). 
By Property 3 of Lemma 3.2 we have, 

sA^{l-s)BCtA^{l-t)B , s<t . (4.8) 
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From (4.7) we obtain, 



(tA0(l-t)i?) -M(s^0(l-s)i?)| \t-s\it,{A)-f,{B)) , 



which in view of (4.8) and (2.1) leads to the proof of the submetric property. 

To prove the measure property, we have to show that for t e [0, 1] (4.5) holds. Observe from Property 1 
of Lemma 3.2 that, 

hiO)=^^iB),hil)^^i{A) , (4.9) 

thus by the monotonicity of h, 

mA,B <tiiB)<ti (A) < Ma,b , 

and for t e [0, 1], 

rriA.B < M (S) < tfi (A) + {I - t) fi{B) < ii (A) < Ma,b- 

The metric property follows from the measure property, (4.8) and (2.1). Finally the inclusion property, follows 
from the assumption B C A, the interpolation property and (4.8). □ □ 

Next we define the measure average in case of simply different sets A, B that do not satisfy the zero- 
measure condition. In view of Remark 3.7, this case is non-generic and in applications can be resolved by a 
small perturbation of the input sets. For completeness we provide a formal construction treating this case. 

An r- offset of a set B with r > 0, is defined as, 

0{B,r) = {p:d{p,B)<r} . 

In case of A, B not satisfying the zero-measure condition, we intersect the set A with an r-offset of B, where 
r is chosen so that the measure of the intersection equals the average of the measures of A, B. 

Definition 4.4. Let A, B be simply different sets that do not satisfy the zero-measure condition. For t £ [0, 1] , 
the measure average of A, B is defined by, 

tA0 (1 - t) B = ci (O {B, rAMt)) n ^) ' (4.10) 

where rA,B{t) is any number in the collection, 

[r■.^^(^0{B,r)f]A)^t^^{A) + {l-t)^l{B)} , (4.11) 

which is not empty, as is proved in Lemma 4-5. For t ^ [0, 1] the measure average is defined as in Definition 
12. 

Note that for any ri,r2 in the collection (4.11), 

ci (O {B, n) fl a) = ci (O {B, r2) fl a) , 

so any r in the collection (4.11) can be used in (4.10). 

In the next lemma we show that the average defined above has the desired properties. 

Lemma 4.5. The measure average of simply different sets in Definition 4-4 well defined and satisfies 
Properties 1-6 in List of Properties 4-i- 

Proof. First we consider the function ifj (r) — fi {O {B, r)f]A). It is easy to observe that, ip (0) = m(-S) and 
that for r large enough ?/; (r) — fJ,{A). To show that ip is continuous, we use a result from [17], guaranteeing that 
for any B CM.^ and X > 0, fi {p : d {p, B) = X) = 0, followed by arguments as in the proof of Lemma 3.8. By 
the continuity of tp, for any t G [0, 1] the collection (4.11) is non-empty. Moreover the set O {B, rA,B{t)) f] A 
is Jordan measurable, consequently in view of (2.2), 

M (^^0 (1 -t)B)^ tfi (A) + {l-t)fi [B) , (4.12) 
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Fig 3. The average tA^{l — t)B of simply different sets for various values of the averaging parameter t. 

implying the measure property. 

Next we observe that the closure, the interpolation and the inclusion properties follow directly from 
Definition 4.4. 

Furthermore the construction in Definition 4.4 yields, 

sA0(l-s)BCtA0(l-t)S , 0<s<i<l. (4.13) 

In all other cases oi s < t, s,t G M., (4.13) follows from (4.8) and the inclusion property of the measure 
average. In view of relations (2.1) and (4.13), the metric property follows from the measure property. 

It remains to prove the submetric property. Let h{x), rriA b and Ma b be defined as in the proof of Lemma 
4.3. From Corollary 3.9, we know that h is continuous anywhere except aX x — ^. From Lemma 3.2, it follows 
that h{0) = IJ,{B), h{l) — IJ,{A), so since h is monotone non-decreasing, h (i) e [fi{B), fi{A)]. We conclude 
that z satisfying, 

mA,B < z < fi{B) or fi{A) < z < Ma,b , (4-14) 

is in the range of h{x). This combined with (4.12) gives a set of relations similar to (4.7), which in view of 
(4.13) and (2.1) leads to the submetric property. □ □ 

Although the List of Properties 4.1 is satisfied by the measure average of two simply different sets satisfying 
or not satisfying the zero-measure condition, the distinction between the two cases is important from the 
geometric point of view. This is so, since the average defined in Definition 4.2 takes into account the geometric 
structure of both sets, while the average in Definition 4.4 is biased towards the smaller set. 

It follows from relations (4.8), (4.13), the interpolation property and the inclusion property, that 

tA0(l -<)BCB, t < , (4.15) 

-BCtA0(l-t)BC A, te [0, 1] , (4.16) 

and 

BC ACiA0(l-i)B, i> 1 . (4.17) 

Relations (4.15)-(4.17), imply that the average tA^{l — t)B "cuts" into the set B for t < 0, interpolates 
between the two sets for t G [0, 1] and expands beyond the set A ioi t > 1. This geometric behavior brings 
the ideas of interpolation and extrapolation into the context of sets. An example of the measure average of 
two simply different sets in ^2 for varying values of the averaging parameter t is given in Figure 3. 

4-2. The measure average of "nested sets" 

We continue the construction with the case oi A,B S Zn such that B C A. We term two such sets nested. 
Formally, the measure average ^ for simply different sets can be straightforwardly applied to any two nested 
sets, preserving all the properties in List of Properties 4.1. However in such a construction, the function g{t) 
in (4.2) dictating the change of the averaging parameter in (4.1) is global. Consequently, when used for 
two nested sets A,B that are not simply different, the reparametrization creates interdependence between 
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different connected components oi A\B, and so may lead to unsatisfactory results from the geometric point 
of view. To reflect the local changes of the geometry of the two sets, we decompose the average of A, B into 
several averages of simply different sets, using B and each of the connected components oi A\B. We denote 
the collection of all connected components of a set D by ^{D). 

Let A, B be nested sets, and assume at first that the number of elements in <t{A \ B) is finite. For any 
C € £(A \ B), we define using the measure average for simply different sets the set, 

Rc,t^t(B[]c)@{l-t)B . (4.18) 

Note that B[_jC and B are simply different sets. 

The results of the averages of simply different sets obtained in (4.18) for all C G £ (A \ B) are merged into 
the average of A and B, taking into account the interpolation and the extrapolation induced by relations 
(4.15)-(4.17). For t > 0, it is logical to take the union of the averages in (4.18), while for i < 0, it is logical 
to remove from B the union of the parts that are "cut" from B by each connected component (see (4.15)), 
which is equivalent to the intersection of the averages in (4.18). We formalize the above procedure in the 
next definition, 

Definition 4.6. Let A, B be nested sets such that the number of elements in €{ A \ B) is finite. The measure 
average of A, B is, 

{ U Ret t^o 

tA^{l-t)B={ ( \ (4.19) 

ci n Ret] t<0 
L \ceciA\B) J 

The operator ci is applied in case t < 0, because regular compact sets are not closed under finite intersec- 
tions. However since P| Rc,t is Jordan measurable, from(2.2) 

Ce€{A\B) 

M n ^c,t = A* ci fl Rc,t 
\c'ec{A\B) J \ \ce(<:iA\B) 

In the study of the properties of the measure average in Definition 4.6, we need the simple observation 
below. 

Lemma 4.7. Let Fi, Fn, Ei, En be sets such that Fi Q Ei,i ^ 1, ...,n, then 

/ n \ / ^ \ ^ 

1- [V}EA\[V}FAQy} {E,\F,) 



2- [i^EA\{^FA(Z{j{E,\F.^ 



Under the stronger assumption that for all i, j £ {1, ...,n} , Fj C E,, 



u ^0 M n ^0 ^ u iE^\F.,) 

i=l / \i=l / i=l 

Proof. Properties 1-2 are immediate to verify. To observe Property 3, let p E ^ jj Ei^ \ ^ p| Fj^ , then 

p E El for some I. Now if p ^ Fi, then p E Ei \ Fi. li p E Fi, then due to the stronger assumption, p E Ei for 
all i. Since there is j such that p ^ Fj, p E Ej \ Fj. DWe can now show that. 

Lemma 4.8. The measure average of nested sets satisfies Properties 1-6 in List of Properties 4.1. 

Proof. The closure property, the interpolation property and the inclusion property follow immediately from 
Definition 4.6 and the corresponding properties of the measure average of simply different sets. Assume 
without loss of generality that s < t, then by relations (4.8) and (4.13), Rc,s ^ Ret, leading to 

s^0(l-s)BCtA0(l-t)B . (4.20) 
12 



To prove the submetric property, we have to consider three cases: (i). < s < t (ii).s < t < and 
(iii).s < < i. Observe that the averages sA^ {1 — s) B, tA^ {1 — t) B in cases (i)-(iii), when written 
in terms of {Rc,s '■ C G €{A\B)} and {Rc,t '■ C G €{A\B)}, correspond to cases 1 - 3 in Lemma 4.7. In 
particular, the assumptions required for case 3 in Lemma 4.7 are satisfied due to relations (4.15)-(4.17). We 
conclude that in view of (4.20) and Lemma 4.7, 

^m(^^0(1-O-B,sA0(1-s)b) =Ai(iA0(l-t)B\sA0(l-s)s) sc:^ I U Rc,t\Rc,s 

\Ceit(A\B) 

Now by the submetric property of the measure average of simply different sets, 

iJi{RcARc.s)< \t-s\^l{C) = \t-s\{^l{A)-^x{B)) , 

C&€(A\B) C£€(A\B) 

yielding the submetric property. 

To prove the measure property we take t £ [0, 1], and observe that for Ci, C2 G (i{A\B), Ci ^ C2 we 
have by (4.16) that Rc^^tClRcz^t = B. Therefore, 




= fi{B)+ Y ^ \ B) , 

Ce£{A\B} 



and by the measure property of the measure average of simply different sets, 

/i(tA0(l-t)s) =A^(B)+ Y tfi{C)=fi{B)+t(p(A)~f,{B)). 

Ce€{A\B) 

Finally the metric property follows from the measure property, by (2.1) and (4.20). □ □ 

The measure average of nested sets satisfies relations (4.15)-(4.17) as well. 

Although the measure average of A with itself is not defined by Definition 4.6, it follows from continuity 
arguments. 

Remark 4.9. If follows from the interpolation property and the metric property, that the sequence of mea- 
sure averages tAi ® (1 — i) B , with Ai = B[J Di and lim /i [Di] — 0, satisfies lim tAi ^ [1 — t) B ~ B . 

Therefore by continuity, we define for any < G M, A G Zn, tA^ {1 — t) A ^ A. 

For the sake of completeness, we consider the case of nested A, B, when A \ B has an infinite number 
of connected components. In this case, since the measure of A is bounded, there is only a finite number of 
connected components of A with measure greater that a preassigned e > 0, and all connected components 
oi A \ B with measure smaller than e are joined into one set, 

UA,B,e = U C . 

CeCiA,B),ti.{C)<e 

The set UA,B,e is treated as a "single component" in Definition 4.6. It is not difficult to show that B [J UA.B,e 
is in 2n- One can verify that all properties in List of Properties 4.1 are preserved in this case. 

4.3. The measure average of general sets 

We are now in position to define the measure average tA@ [1 — t) B of two general sets A,B £ The 
average is decomposed into two averages of nested sets, 

i?i,t = iA0(l-i)ci(Af|B) , (4.21) 

and 

-R2,t = (l-i)S0ici(Af|B) . (4.22) 

The two averages are merged preserving the geometry of the interpolation and the extrapolation of sets. 
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Definition 4.10. Let A, B G Zn, the measure average of A, B with < G M is, 

' Ri,t[jR2,t te [0,1] 

a0(l-i)B=i ci(i?i,A^n5)U^2,t t>l (4.23) 

{l-t)A^tB t<0 
with the sets Ri.t,R2.t defined in (4-.21)-(4-22). 

Notice that ior t > 1, Af] B is removed from Ri^t, so that the "cutting" from Af]B hy the extrapolation 
in i?2,t will affect the resulting average. 

Remark 4.11. It follows from Definition 4. 10 that for any t e R, A, B e a„, tA^{l-t)B = (1 - t)B^tA. 

Theorem 4.12. The measure average of any two sets A, i? e 3™ satisfies Properties 1-6 in List of Properties 
11. 

Proof. The closure, the interpolation and the inclusion properties follow from the similar properties of the 
measure average of nested sets. Next, we observe that from the inclusion property in the nested case, we 
have for i e [0, 1], 

i?i,tn^2,t = ci(Af|B] , (4.24) 

so, 

^i(tA^{l-t)B)=^i [Ri^t] + M (i?2,t) - M (a fl . (4.25) 
Using (4.25) and the measure property in the nested case we obtain that for t e [0, 1], 

^^(tA@{l-t)B)^t^Ji{A) + {l-t)^i(A[]B)+{l-t)^i{B) + t^l{A[]B)-^i(A[]B) , 

which yields the measure property. The metric property is proved using (4.24) and the metric property of 
the measure average in the nested case. 

To prove the submetric property, we observe that for s <t, i?i.s C Rij. and R2.t ^ R2.s- Consequently, 

dt,[sA^{l~s)B,tA^{l-t)B^ <df,iRi,„Ri,t) + d^,iR2,s,R2^t) , (4.26) 

leading to the submetric property in view of the submetric property of the measure average in the nested 
case. □ 

□ 

As an immediate consequence of the metric property of the measure average, we obtain a result about 
the approximation of Holder continuous SVFs by the piecewise interpolant based on the measure average. 
We term a SVF F H61der-z^ continuous if, 

d^{F{ti),F{t2))<C\ti-t2\'' , (4.27) 

where C is a constant (termed the Holder constant of F) depending on F but not on ^1,^2 and e (0, 1]. A 
Holder- 1 continuous SVF is also termed Lipschitz continuous. 

Corollary 4.13. Let F : [0, 1] — ?> 5„ be H6lder-v continuous. We define P/y {x) ■ [0, 1] to be the 

piecewise interpolant, 

Pn (x) = - «) ^ (*^) ((* + 1) - ^) -P" ((« + 1) > 2; e [th, {i + l)h], h = 1/N, i = 0, ...N + I . 

(4.28) 

Then for any x G [0, 1], 

d^ {F{x),PN{x))<Ch'' , 

where C is the Holder constant of F. 

In the following sections we study subdivision methods based on the measure average of sets. 
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5. Spline subdivision schemes adapted to sets with the measure average 

In this section, we use the measure average in the adaptation to sets of sphne subdivision schemes (see.e.g [14], 
Section 3.1) . First, the refinement rule is expressed by repeated binary averages using the Lane-Riesenfeld 
algorithm [24]. Binary averages of numbers are then replaced by the measure averages of sets. Our approach 
is similar to [11], where spline subdivision schemes are adapted to sets using the metric average. 
An m-degree spline subdivision scheme (m > 1) in the real-valued setting refines the numbers, 

{ft : I e Z} C M, 

according to the refinement rule, 

/'^'=E4'^]+.-2/i^ ,fe = 0,l,2,..., (5.1) 



where [•] is the ceiling function, a^™^ — ^ ^ ^ Z^™ for I = 0,1,..., rn + 1, and al™-* = for / G 

Z\{0,l,...,m + 1}. 

The case m = 2 is the Chaikin subdivision scheme with the refinement rule, 

fk+l _ f fc I 1 i-fe 



•'21+1 ~ A Ji A -li+l ■ 



-f'+'- 

4-'' ^ 4'' 

Chaikin scheme is the simplest which generates limits, and is widely used. 

At each level k, the piecewise linear interpolant to the data (2~'^i, f^) , fc e Z is defined on M by. 



/fe(x) = A(a:)/f + (l-A(x))/ti , 12-" <x<{i + l)2-\ (5.2) 

where A (x) = {i + 1) — x2^ . The sequence {/^ (x)} converges uniformly to a continuous function f°° (x), 
which is the spline of degree m with the control points {i, fi) ,i G Z. 

The Lane-Riesenfeld algorithm evaluates (5.1), by first doubling the values at level k and then performing 
m steps of repeated averages. We apply the above procedure to sets, by replacing averages of numbers with 
the measure average of sets. We combine the doubling step with one averaging step. So first the sequence of 
sets at level fc, {F/^' : i e Z} is refined using the measure average, 

F^r^'^Ft, ^tti" = ^^?0j^ti, ^eZ. (5.3) 



Then for 1 < j < m — 1, the sequence is replaced by the measure averages of pairs of 

consecutive sets. 



pk + l.j 



\Fr'^'-' ®\Ft:^^^-\ zeZ. (5.4) 



2 

Finally, the refined sets at level fc -I- 1 are given by, 

F^^' = , ^eZ. (5.5) 

For the analysis of the above family of subdivision schemes we use, similarly to the real-valued case, the 
piecewise interpolant Fk (x) in terms of the measure average through the sets al the k-th level, 

Fk{x)^Xix)Ft^{l-\ix))Ft:,, , z2-'=<a;<(z + l)2-\ (5.6) 

where A (a;) — {i + 1) — x2^ . 

The following results on convergence and approximation order of spline subdivision schemes based on the 
measure average are analogous to results on spline subdivision schemes bases on the metric average in [10]. 
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Theorem 5.1. The sequence of set-valued functions {Fk (a;)}^,^^^ converges uniformly to a continuous set- 
valued function F°° (x) : M — >■ £„, which is Lipschitz continuous relative to the symmetric difference metric 
with the Lipschitz constant L — supd^ {F^ , F^_^_i^ . 

i 

Theorem 5.2. Let G : M — >■ £„ be Lipschitz continuous, and let the initial sets be given by F.^ = G {S + ih) G 
3n; i G Z with 5 e [0, h) and h > 0. Then, 

d^{G{x),Fk{x))<Ch , (5.7) 

where F^ (x) is given by (5.6) and C is a constant depending on the degree of the scheme. 

Corollary 5.3. Under the assumptions of Theorem 5.2, the distance between the original set-valued function 
G {x) and the limit set-valued function F°° (x) is bounded by, 

m&xd. {F°° {x),G{x)) < Ch . 

X 

The properties of the metric average relative to the Hausdorff metric, used in the proofs of results analogues 
to the above results, are also possessed by the measure average relative to the symmetric difference metric. 
Therefore the proofs of Theorems 5.1,5.2 and Corollary 5.3 are similar to the proofs of Theorems 4.3, 
4.4 and Corollary 4.5 in [10] respectively, and are omitted here. Theorem 5.2 and Corollary 5.3 can be 
straightforwardly extended to H61der-z^ SVFs to obtain approximation order O (ft."). 

Next we use the measure property and the inclusion property, which are specific to the measure average, 
to derive further properties of spline subdivision schemes based on the measure average. As a consequence 
of the measure property we have. 

Corollary 5.4. Let S be an averaging rule defined for a sequence {fi}^^^ C R by, 

S ({/iLez) = Y^a^f, , 

i 

with Oi > 0, ^ fli — 1. Let S* be an adaptation of S to sets by representing S as a sequence of repeated 

i 

binary averages of numbers with non-negative averaging parameters, and replacing averages of numbers by 
the measure averages of sets. Then for any sequence of sets {-Ffjjg^ 

^^{s*{{F,}^^^))^s{{^i{F,)}^^^) . 

Note that the result of Corollary 5.4 is independent of the specific representation of the averaging rule S 
by repeated binary averages. 

It follows from Theorem 5.1 and the definition of the metric (•, •) that, 

tJi{F^{x)) = lim ^l[Fk{x)) . (5.8) 
Corollary 5.4 together with (5.8) leads to. 

Corollary 5.5. . Let {^f'l^gg C 5n. Let F°° be the limit SVF of the set-valued spline subdivision scheme 
applied to {F^}^^^ and let f°° be the limit function of the real-valued spline subdivision scheme applied to 
{l^{F")}.^z- Then, 

/.(F-(:r)) = /-(x) . 

Next we state several results concerning set-valued spline subdivision schemes applied to monotone data. 
It is well knows that real- valued spline subdivision schemes are monotonicity preserving. Due to the inclusion 
property of the measure average and relations (5.3)-(5.5), spline subdivision schemes adapted to sets with 
the measure average are also monotonicity preserving. A sequence of sets {Fijjgx is termed monotone 
non-decreasing (non-increasing) if Fi C Fi^i {Fi I) -Fi+i). With a similar definition for a non-decreasing 
(non-increasing) set valued function we obtain, 
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Corollary 5.6. Let {Fi}.^^^. monotone non-decreasing (non-increasing). Let F°° be the limit 

SVF of the set-valued spline subdivision scheme applied to {Fi}^^^^- Then F°° is monotone non- decreasing 
(non-increasing) . 

The notion of the speed of a curve in a metric space {see e.g. [7], Chapter 2), can be used as an indication 
of the "smoothness" of a set- valued function. For a real- valued / the speed at a point x is, 

, , ,. \fix)-f{x + e)\ 
Vf (x) = lim 



whenever the limit exists. For differentiable /, Vf is the absolute value of the derivative of /. We define the 
velocity of a SVF F, 



vp (x) = lim 



d^ (F {x) ,Fix + e)) 



By combining relation (2.1) with Corollaries 5.5 and 5.6 we arrive at, 

Corollary 5.7. Let {-Ffj^g^ monotone non- decreasing (non-increasing). Let F°° and f°° he defined 

as in Corollary 5.5, then 

Vp'^ (x) = (x) . 

Under the assumptions of Corollary 5.7, we have that for the spline subdivision of degree m > 2, vpoa is 
continuous and has continuous derivatives up to order m — 2. 

In the next section we adapt to SVFs the 4-point interpolatory subdivision scheme using the measure 
average with both positive and negative averaging parameters. 

6. The 4-point subdivision scheme adapted to sets with the measure average 

In the real-valued setting, the 4-point subdivision scheme is defined by the following refinement rule 

&'^ft, (6.1) 

and 

&+\ = (/ti + fU + (1/2 + uj) (/f + /ti) , (6.2) 

repeatedly applied to refine the values {//^ : i e Z} c M for fc = 0, 1, 2.... Here w is a fixed tension parameter. 
Usually w is chosen to be 1/16, since this value yields the highest approximation order [15] and the maximal 
Holder exponent of the first derivative of the limit function [8]. The rule can also be applied to refine a 
sequence of points in R". 

One can see that the coefficients in (6.2) sum to one, so it is a weighted average of the four values 
fi+2- Rewriting the insertion rule (6.2) in terms of binary weighted averages as, 

&+\ = ^((-2u;)/^l + (l + 2u;)/f) (6.3) 

+ ^-{{-2w)f^+, + {l + 2w)ft+,) , 

we adapt the refinement rule defined by (6.1) and (6.3) to sets by replacing binary averages of points by the 
measure averages of sets. 

The refinement rule for the sets {F^ : i g} C 5n is. 



F^+' = Ft , (6.4) 



and 



^^\ = ^^^0^^> (6.5) 
with 

= {-2w) Ft, 0(1 + 2w) Ft , (6.6) 
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and 

H^+\ = (-2^) Ft+2 0(1 + 2w) Ft+, . (6.7) 

Note that the subdivision with the refinement rule (6.4)-(6.5) is interpolatory. 
First we study the convergence of the scheme. 

Lemma 6.1. Let {Pj; : i e Z} 



ne, 



dk = supdf,{F^^,FlXi) , 



then 



dk < • (6.8) 

Proof. By the interpolation property and the submetric property of the measure average (see List of Prop- 
erties 4.1), 

{E2^+\ , Ft) < ^y^d^ {Fl,, Ft) < 2wdk , (6.9) 

and, 

d^. {H^^^+\,Ft+i) < 2wd^ {Ft+i^F^^^) < 2wdk . (6.10) 

Therefore, 

d, {E^2t+i,H^2t+\) < {E'2t+\,Ft) + {Ft,FlX^ + d, {H^+\, Ft,,) < (1 + Aw) du , (6.11) 
from which we obtain by using (6.5) and the metric property, 

d, {F2%\ , E^+,\) < Q + 2z«) 4 . (6.12) 

Thus by the triangle inequality we get from (6.9) and (6.12), 

dAF2T,F,%\)<(^l+4w^ dk. 

Similarly one gets, 

dAF2tlF2%\)<(^l+Aw^ dk. 

Therefore 

dk+i - Q + ^"^^ ' 

and (6.8) holds. □ □ 

We conclude from Lemma 6.1, 
Lemma 6.2. Let {Fk {x)}|^^^_^_ be the sequence of piecewise interpolants defined as in (5.6), then 

d^ {Fk (x) , Fk+i (x)) < C Q + 4u;^ , (6.13) 

with C = do (1 + 4u;). 

Proof. For i2-^ < x < {i + \) 2'^ , using Ft = F^t^ , 



d^ [Fk (x) , Fk+i (x)) < d^ {Fk (x) , Ft) + d^ {Ft+\Fk+i {x)) < ^k + dk+i 
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and in view of (6.8), 

d^, {Fk {x) , Fk+i (a;)) < ^do ( ^ + 4w ) +do[^+Aw 



leading to the claim of the lemma in this case. A similar argument for (i + 2 < x < {i + 1)2 with 
Fj^ replaced by Ff^^^ completes the proof. □ □ 

Theorem 6.3. The sequence of set-valued functions {Fk {x)}k^x+ converges uniformly to a continuous set- 
valued function F°° (x) : M £„ whenever w < 

Proof. By definition the functions {Fk {x)}k^z+ continuous. By the triangle inequality, 

k+M-l 

{Fk (x) , Fk+M (x)) < ^ {F, (x) , ^^,+1 {x)) . 



i—k 



From Lemma 6.2 and by the assumption w < 



8 ' 

df,{Fk{x),Fk+M{x)) < do{l + 4w)[-+4w\ < do {1 + Aw) i- + 4w j • (6-14) 



i—k 



We observe from (6.14) that for w < |, {Fk{x}}i.^^^ is a Cauchy sequence in the metric space {3n,d^}, 
and consequently it is also a Cauchy sequence in the metric space of Lebesgue measurable sets with the 
metric (•, •), {£„, c?^}. Since d^} is a complete metric space, the sequence {Fk {x}}k^z converges to 
F°° (x) e £„■ The convergence is uniform in x due to (6.14), consequently F°° (x) is continuous. □ □ 

Next we derive results concerning approximation of Holder continuous SVFs by the 4-point subdivision 
scheme. 

Theorem 6.4. Let G :M. ^ J„ be v-Holder continuous, and let the initial sets be given by F^ = G {6 -\- ih) 
i e Z with S e [0, h) and h > 0. Then for Fk (x) given by (5.6), 

d^{G{x),Fk{x))<Ch- , (6.15) 
where C = + i + 1/2-^4'tu ) ^ '^^'^ ^ Holder constant of G. 

Proof. Without loss of generality assume that 5 = 0. Let x be such that ih < x < (i + 5) /i. From the 
triangle inequality, 

d^ (G (x) , Fk (x)) < d^ (G (x) , F,) + d^ {F,, Fo (x)) + d^ (Fq (x) , Fk (x)) . 

Using the Holder continuity of G, the metric property of the measure average and (6.14) we get. 



df,iG{x),Fk{x))^H(^^h^ +ldo + l±l^do. (6.16) 



Now c?o < Hh'^ with H the Holder constant of G, and therefore (6.16) implies (6.15). □ □ 

Corollary 6.5. Under the assumptions of Theorem 6.4, the distance between the original set-valued function 
G (x) and the limit set-valued function F°° (x) is bounded by, 

maxd^ {F°° (x) , G (x)) < Ch" , 

with G as in Theorem 6.4 ■ 

Since in extrapolation the measure property does not hold, there is no result for the 4-point scheme similar 
to that in Corollary 5.5 for spline subdivision schemes. 
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7. Computational examples 

The algorithms for the computation of the measure average along with implementation details and the 
numerical evaluation of our methods applied to the reconstruction of 3D objects from cross-sections are given 
in [23]. Here we provide several computational examples illustrating the theory presented in the previous 
sections. 

In our examples we use the adaptations to sets with the measure average of the piecewise linear interpo- 
lation (c./. Section 4), Chaikin subdivision scheme (c./. Section 5) and the 4-point subdivision scheme (c./. 
Section 6). Chaikin scheme is implemented using the Lane-Riesenfeld algorithm as described in Section 5. 

While the subdivision schemes discussed so far are defined over the whole real axis, computations require 
to deal with a finite number of sets, ant therefore to consider boundary rules for these schemes. Boundary 
rules for the real- valued subdivision schemes are discussed in ([29], Chapter 32). 

Assume that at the k-th level of the subdivision, we have the sets F^,..., F,j , assigned to equidistant points 
Xq, a;^. For Chaikin scheme, the refinement rules at the boundaries are, 

pfe+l 1^ T^k /T\ I pk 

^1 ^ 2 tP 2 1 ' 

pfc+1 _ 1 p/c rr\ 1 pfe p'^+i pfc 

^2n-2 ^ 2 n-1 2 " ' ^ 2n-l ~ ^11 1 

and FJ', i — 2, ...,2n — 3 are computed using relations (5.3)-(5.5). With these boundary rules the limit of 
the Chaikin scheme interpolates the sets at the endpoints. For the 4-point scheme, the modified refinement 
rules at the boundaries are, 

pfc+i 1 pfc rr\ 1 pfe. — 1 pfe n:^ i pfc 

while F^+^, i = 0, 2, 3, 2n - 3, 2n - 2, 2n are obtained using relations (6.4)-(6.5). In both schemes, the 
refined sets at level k+\ are assigned to points, which are obtained from Xq, x^ by applying rules analogous 
to those applied to the sets. 

In the first example we apply the three aforementioned methods to a collection of eight sets Fq..., i^y with 
varying topologies, located on equidistant parallel planes. The results are visualized in Figures 4-6. Note 
that in all examples only boundaries of the sets are depicted. It can be observed from Figure 4, that the 
piecewise interpolation indeed passes through the original sets, but the transitions between pairs of original 
consecutive sets are noticeable. Figure 5 demonstrates the smoothing effect of Chaikin subdivision, however 
the limit function does not pass through the original sets. Finally, the limit of the 4-point subdivision scheme 
in Figure 6 has a certain smoothing effect, and yet it passes through the original sets. The above geometric 
behavior of the three set-valued methods is in analogy with the well known features of the corresponding 
real- valued methods. 

In Figure 7 we plot {x)), where {x), i = 1,2,3, are the SVFs in Figures 4-6 respectively. 

In accordance with Corollaries 5.4 and 5.5, /i (F'^I (x)) and /i (F^^I (x)) are equal to the piecewise linear 
intcrpolant and the limit of Chaikin subdivision applied to fi{Fi), i = 0, ...,7. While a result similar to 
Corollary 5.5 does not hold for the 4-point scheme, /x (F3 (x j) appears as a smooth curve interpolating the 
data ^(Fq) , ...,^(F7). 

Next we consider the piecewise interpolation and Chaikin subdivision applied to monotone data Fg, F7, 
Fi 3 Fi_|_i, i = 0, 6. In view of the inclusion property of the measure average and Corollary 5.6 the resulting 

SVFs are monotone in the set-valued sense, as illustrated in Figure 8. In Figure 9 we plot fi ^F^*! {x)j , where 

fW (x), i = 1,2, are the SVFs obtained using piecewise interpolation and Chaikin subdivision respectively. 

By Corollary 5.5, /i ^F^*! {x)j , i = 1,2 are the same as the functions obtained by the application of piecewise 

linear interpolation and Chaikin subdivision respectively to the initial data fi (-^i) i * = Oi 7. In particular 
/i ^F^^l (x)^ has a continuous derivative, which in view of Corollary 5.7 illustrates the continuity of the metric 
velocity of F^^l (x). 



-^0 ^ -^0 ' 
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a b c 




Fig 4. Piecewise interpolation of sets: a. the initial sets; b. 15 sets introduced between each pair of consecutive sets using the 
measure average; c. visualization of the resulting SVF as a 3D object. 



a b c 




Fig 5. Refinement of sets with Chaikin subdivision scheme: a. the initial sets; b,c,d,e. the refined sets after 1,2,3,4 subdivision 
steps respectively; f. visualization of the resulting SVF as a 3D object. 



8. Extensions 

We extend some of the results obtained in metric spaces of sets in this work and in [10] to general metric 
spaces endowed with a binary average satisfying certain properties. Let {X,dx} be a metric space, and let 
ffl be an average on elements of X defined for non-negative averaging parameters {S : [0,1] x X x X ^ X). 
Assume that the average ffl satisfies the interpolation property and the metric property in List of Properties 
4.1 with -j;„, df^ and replaced by X,dx and ffl respectively. Then: 

1. A piecewise interpolant based on ffl can be defined as in (4.28), leading to an approximation result 
analogous to Corollary 4.13. 

2. The spline subdivision schemes can be defined using the Lane-Riesenfeld algorithm. Such an adaptation 
leads to convergence and approximation results analogous to Theorems 5.1,5.2 and Corollary 5.3, with 
the limit of the subdivision in the metric completion of {X, dx}- 

3. Under the additional assumptions that the average ffl is defined also for averaging parameters outside 
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d e f 




Fig 6. Refinement of sets with the 4-point subdivision scheme: a. the initial sets; b,c,d,e. the refined sets after 1,2,3,4 subdivision 
steps respectively; f. visualization of the resulting SVF as a 3D object. 




Fig 7. (_fW (x)), i = 1,2,3, for the SVFs in Figures 4, 5, 6, depicted by the dash-dotted, the dotted and the solid lines 
respectively. The measures of the initial sets are denoted by * . 



[0, 1] and satisfies the submetric property in List of Properties 4.1,ttie 4-point subdivision scheme can 
be adapted to the elements of X, as in relations (6.4)-(6.5). Convergence and approximation results 
analogous to Theorems 6.3, 6.4 and Corollary 6.5 are obtained in a similar way. 
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